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A possible equivalence of scalar dark matter, the infiaton, and modified gravity 
is analyzed. After a conformal mapping, the dependence of the effective Lagrangian 
on the curvature is not only singular but also bifurcates into several almost Ein- 
steinian spaces, distinguished only by a different effective gravitational strength and 
cosmological constant. A swallow tail catastrophe in the bifurcation set indicates the 
possibility for the coexistence of different Einsteinian domains in our Universe. This 
'triple unification' may shed new light on the nature and large scale distribution not 
only of dark matter but also on 'dark energy', regarded as an effective cosmological 
constant, and inflation. 



1. Introduction 

The origin of the principal constituents of the Universe in the form of dark matter 
(DM) and dark energy (DE) (or rather dark tension due to its negative pressure), 
remains a major puzzle of modern cosmology and particle physics. Ad hoc proposals 
like modifying Newton's law of gravity as in Milgrom's MOND (Modified Newtonian 
Dynamics) are difficult to reconcile with relativity, or need a phase coupling to 
a complex scalar. 3 The particle physics' view is to leave gravity untouched but 
postulate WIMPs (Weakly Interacting Massive Particles) like axions, dilatons or 
neutralinos as dark matter candidates. 66 

The dominant non- visible "dark" fraction of the total energy of the Universe is 
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known to exist from its gravitational effects. Since the dark matter part is distributed 
over rather large distances, its interaction, including a possible self- interaction, 60 
must be weak. The main candidates for such weakly interacting particles are the 
(universal) axion or the lightest supersymmetric particle, as the neutralino in most 
models, cf. Ref. 66 

Recently, an observed excess of diffuse gamma rays has been attributed 8 to the 
annihilation of DM in our Galaxy. The flux of the presumed neutralino annihilation 
allows a reconstruction of the distribution of DM in our Galaxy. Most probable is a 
pseudo-isothermal profile but with a substructure of two doughnut-shape rings in 
the galactic plane. It is believed that these transient substructures have their origin 
in the hierarchical clustering of DM into galaxies. However, there are reservations 
concerning the internal consistency in the interpretation of the observations: Given 
the same normalization for the cross section of the DM particle, it appears unlikely 1 
that DM annihilation is the main source of the extragalactic gamma-ray background 
(EGB). Moreover, then the large accompanied antiproton flux needs to meet several 
constraints. 5 

On the other hand, heterotic string theory provides a very light universal axion 
which may avoid 12 the strong CP problem in quantum chromodynamics (QCD). 50 
Given the existence of such almost massless (pseudo-)scalars, it has been speculated 
that a coherent non-topological soliton (NTS) type solution of a nonlinear Klein- 
Gordon equation can account for the observed halo structure, simulating a Bose- 
Einstein condensate of astronomical size; 16 ' 17,34 cf. Ref. 53 In particular, a $ 6 toy 
model 43 yields exact Emden type solutions in flat spacetime, including a flattening 37 
of halos with ellipticity e < I as observed via microlensing. 14,54 

Both views are to some extend physically equivalent, i.e., scalar dark matter 
minimally coupled to Einstein's general relativity (GR) is equivalent 4,27 ' 49 ' 55 ' 59 to a 
modified gravity in the relativistic framework of higher-order curvature Lagrangians. 
Such effective Lagrangians may also arise from the low-energy limit of (super-) 
strings. The model proposed by Carroll et al. also uses such an equivalence of a 
nonlinear higher-order curvature Lagrangian, in order to explain the present cos- 
mic acceleration, but takes resort to an 1/R type curvature Lagrangian, which is 
unbounded for weak gravitational fields. 

Quite generally, the dependence of the effective Lagrangian L on the scalar 
curvature R is singular and cannot always be represented analytically in the (i?, L) 
plane. More stunning is our recent finding 55 that, although nonlinear, a L = L{R) 
type Lagrangian bifurcates in several branches of almost linear Lagrangians 

Lc^-^-(R-2A cS ) (I) 

ZK e ff 

distinguished only by the effective gravitational constant « c ff and a cosmological 
constant A c ff. 

This indicates a unifying picture of dark matter, 'dark energy' 53 with an equation 
of state parameter w DE = p/p ~ —1, and modified gravity which may account, on 
different scales, for inflation, 37 ' 40 dark matter halos 36,43 of galaxies or even dark 
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matter condensations, the so-called boson stars 23 ' 24,41 ' 58 as candidates of MACHOs 
(Massive Compact Halo Objects). The 'landscape' view 25 of (super-)strings suggests 
also such a triple unification based on coherent, possibly oscillating scalar field with 
a mass larger than the Hubble scale at the present epoch, i.e. m > = 10~ 23 eV. 

2. Gravitationally coupled scalar fields 

Conventionally, dark matter and inflation can be modeled by a real (or complex) 
scalar field cj> with self-interaction /7(|0| 2 ) minimally coupled to gravity with the 
Lagrangian density 



C 



DM 



{R + «[> - 2t/(|0| 2 )] } , (2) 



2 re 

where re — 8nG is the gravitational constant in natural units, g the determinant of 
the metric g M „, and R the scalar curvature of Riemannian spacetime with Tolman's 
sign conventions. 63 A constant potential Uo = A/re would simulate the cosmological 
constant A. 

In a wide range of inflationary models, the underlying dynamics is simply that 
of a single scalar field — the inflaton — rolling in some underlying potential. This 
scenario invented by Linde 26 is generically referred to as chaotic inflation due to its 
choice of initial conditions. Many superficially more complicated models can also 
be rewritten in this framework. 



3. General metric of a flat inflationary universe 

A spatially flat {k = 0) Friedman-Robertson- Walker (FRW) Universe with metric 

ds 2 = dt 2 - a 2 {t) [dr 2 + r 2 (d9 2 + sin 2 6dip 2 )} (3) 

is nowadays favored by observations. 61 Its temporal evolution of the generic model 
l|2|) is determined by the two autonomous first order equations 



H 



-3ff 2 (l 



wde) = k.U{4>) - 



:m 2 



V(H,< 



re 



(4) 
(5) 



where H = d(t)/a(t) is the Hubble expansion rate. 

The function V(H, </>) will turn to be the "height function" in Morse theory. 46 
Observe that V < in order to avoid scalar ghosts. For the FRW metric, the 
Lagrangian density j2|) reduces to 



C = -■ 



-d 2 a + 



j, 2 - u{4>) 



(6) 



Since the shift function is normalized to one for the FRW metric, the canonical 
momenta are given by P = dC/dd = — QHa 2 /n and it = dC/d<j) = a 3 cp. The 
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Hamiltonian or "energy function" is given by 

E 



U 2 + -V{H,< 

I K 







(7) 



and, due to vanishes for all solutions, which is a familiar constraint in GR. 

Using the Hubble expansion rate H := a(t)/a(t) as the new "inverse time 
coordinate 

dH 



t = t(H) = f - 

J K 



U-3H 2 



we were able to find the general solution 



r>(> 



a = a(H) = ao cxp ( / — 



11,111 \ 



U-3H 2 



ds 2 = 



dH 2 



kU - 3H 2 



qq exp 2 



HdH 

kU - 3H 2 



[dr 2 ^ 



r 2 (d6 2 + sin 2 6d(p 2 )] , 



dH 



V'SH 2 - kU 



(8) 
(9) 

(10) 
(11) 



where U — U(H) :— U(4>{t{H))) is the reparametrized inflationary potential. 

The singular case U = 3H 2 /n with wde = — 1 leads in ((4]) to the de Sitter 
inflation with exponential expansion a(t) — ao exp(y / A/3t), which for a later time 
becomes physically unrealistic, since the inflationary expansion eventually needs to 
merge into the usual Friedman cosmos. Therefore, in explicit models we use the 
ansatz 



U(H) = -H' 2 



9(H) 



(12) 



for the potential, where g{H) = V (H, 4>{t(h))) is a nonzero function which should 
provide the graceful exit from the inflationary phase to the Friedman cosmos. In 
order to have a positive acceleration during the inflationary phase, but on the other 
hand a real scalar field in (fTT|) . its allowed range is — H 2 < g < 0. This if-formalism 
facilitates considerably the reconstruction of the inflaton potential. 36 ' 40,56 In second 
order this is based on an Abel equation 13 for the primordial density e = —g/H 2 . 



3.1. Classification by catastrophe theory 

A general classification of all allowed inflationary potentials and scenarios has al- 
ready been achieved by Kusmartsev et al. 22 via the application of catastrophe theory 
to the Hamilton-Jacobi type equations (J4j) and |5]) regarded as an autonomous non- 
linear system. 

In phase space, the equilibrium states are given by the constraint {H, ((>} = 0. 
The critical or equilibrium points, respectively, of this system are determined by 
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V(H c ,(j> c ) = 0. This constraint is globally fulfilled by kU(4>) = iH\, where the 
Hubble expansion rate is constant, i.e. H\ =: \/A/3 corresponding to wde = — 1- 
For = and A 7^ 0, we recover the de Sitter inflation. 
In general, the Jacobi matrix 

_ / -6H kU' 

\±6H(-2kV)- 1 / 2 tkU'(-2kV)-^ 2 

of the system (0J and depends crucially on U' = dU/dtfi. Since its determinant 
vanishes, i.e. det J = 0, the system is always degenerate. For the analysis of stability, 
it suffices therefore to consider only Q and to reconstruct </> later. Since H and <j> 
are independent variables, we can introduce the non-Morse superpotential W(H, cj>), 
defined via V := —dW/dH, and analyze the system with the aid of catastrophe 
theory. From ([¥]) we obtain 

W(H, dp) = H 3 - nU{(j})H + C{(j>) , (14) 

where C is an arbitrary function of (p. The function W is already in canonical 
form in -H-space, and belongs to a Whitney surface 2,21 or to the Arnold singularity 
class Ai. The corresponding Whitney surface has only one control parameter, the 
potential U. Thus, an evolution of critical points is determined via the values of the 
potential U . Let us analyze the types of critical points at different fixed values of 
the control parameter U . 

If U c := U(4> c ) < 0, the equation has no stable critical point due to the shape of 
the Whitney surface. However if U c > 0, there are two critical points: a stable one 
at H c = y/Uc/3 and an unstable one at H c = — ^U c /3- 

The value U c = is the bifurcation point. Provided this is also an extremal of 
V, we necessarily have dV/dH | c = —6H C = 0, dV/d(f> | c = kII' c = 0, and <p c — 0. 
Thus, also the critical points of the Klein-Gordon equation for the scalar field are 
involved. Hence, the Hubble parameter has to vanish and </> c is a double zero of the 
potential U. The Hessian of (@| takes the form 

Hess ^ = ("o 6 ^)- (15) 

The sub-determinants of the Hessian are Ao = d 2 V/dH 2 = —6 < and Ai = 
detHess(V) = {d 2 V/dH 2 ) (d 2 V/d<j) 2 )-{d 2 V/{dHd<j))) 2 = -6kU". For a maximum 
of the potential U, i.e. U' = and U" < 0, the function V possesses a maximum; a 
minimum of the potential U, however, corresponds to a saddle point of V. 

3.2. Reheating 

Since U{4> c ) can be associated with the latent heat of the Universe in this epoch, we 
could demonstrate: 22 

The critical points of the non-Morse potential W{H, 4>) determine the evolution 
in the inflationary phase. Along the minima and maxima H c = ±y/U(<fi c )/3, the in- 
fiaton moves from the slow-roll to the hot regime. The saddle points ofW, i.e. more 
precisely, the minima of V , determine the onset of reheating. 
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3.3. Scale-invariant spectrum as a limiting point 

It is now possible to relate the "slow-roll" condition, for the velocity of the infla- 
tionary phase, to the critical points resulting from catastrophe theory. For inflation 
(with a > 0) the two "slow-roll parameters" are, in first order approximation, 

6 = -g/H 2 and T) = -dg/d(H 2 ) , (16) 

where g = V(H, <fi) is the "graceful exit function" . In this reduced dynamics, 
they are effectively determined by the first and second derivatives of the re- 
duced non-Morse function W(H), i.e., more precisely by e = (1 / H 2 ) (dW/dH) 
and r\ = (l/2H)d 2 W/(dH) 2 . They will also determine the density fluctuations of 
the early Universe. 

The deviation A := (n — l)/2 from the scale-invariant Harrison-Zeldovich spec- 
trum with n = 1 leads in the first order slow-roll approximation to the differential 
equation 

with 

e = -A + AH 2 , V = AH 2 - AH 4 (18) 

as solutions for the density and the graceful exit function, respectively. 
Then the non-Morse potential turns out to be 

W(H, <f>) = - I VdH = -^-H 3 + -H 5 + C{4>) (19) 
J 3 5 

which together with (|14[) provides us with the reparametrized potential 

«£(» = (i + jW--/? 4 , ( 2 °) 

which for A ~ is positive in the leading order. 

In the next order slow-roll approximation arises a nonlinear Abel equation 13 
with a continuous spectrum for n < 1 and a discrete 'blue' one for n > 1, see Ref. 57 
for more details. This division of the spectrum persists in the next to second order 
slow-roll approximation. 36 More important is our finding that in both higher order 
approximations the Harrison-Zeldovich spectrum with n = 1 stays a limiting point, 
which agrees quite well with recent constraints from the observations of WMAP. 20 



4. Higher order curvature Lagrangians via field redefinition 

Nonlinear modifications of the Einstein-Hilbert action are of interest, among oth- 
ers, for the following reasons: First, some quadratic models can be renormalized 
when quantized, cf. Ref. 38 Second, specific nonlinear Lagrangians have the prop- 
erty that the field equations for the metric remain second order as in GR; these are 
the so-called Lovelock actions which arise from dimensional reduction of the Eu- 
ler characteristics, cf. Ref. 47 In Yang's theory of gravity, this topological invariant 
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induces an effective cosmological constant for instanton solutions residing in Ein- 
stein spaces. 33,38 At times, 64 this renormalizable model is referred to as Yang-Mielke 
theory of gravity. 

As an instructive example of higher-derivative theories of gravity, 52 we consider 
the Lagrangian density 

1 „ 1 



C = L(R)y/\ g | = I ^R + ^R 2 j V\gl (21) 

where j3 is a dimensionless coupling constant. Through the conformal change 30 

dL _ ( 2k 
: dR~\ + W 

of the metric, this Lagrangian can be mapped to the usual Hilbert-Einstein La- 
grangian with a particular self-interacting scalar field. Indeed, Starobinsky 62 con- 
sidered earlier such models for inflation. 

The scalar field, the inflaton, will arise via 

/ 3 

4> = J— In ft (23) 

V ZK 

from the nonlinear parts of a higher-order Lagrangian L = L(R) in the scalar cur- 
vature R. Recently, such modified gravity models are considered as alternatives 7 ' 48 
for dark matter or even dark energy. 



la/3 -> g a p = tyap with tt = 2k— = ( 1 + —^R ) , (22) 



4.1. Reparametrized Lagrangian 



Instead of studying the resulting complicated nonlinear field equations of higher- 
order curvature Lagrangians, we follow the equivalence proof of Ref. for the con 
formal frame ([22]) . Then, our inflaton Lagrangian density (J2]) acquires the form 

1 " : 2Kfi 2 u(n) 



c = 



2k 



RQ 



(24) 



where U(Q) :— U((f>(fl)) = U ^-\/3/2rtlnf2j is the reparametrized potential. Thus 
in our approach, the inflaton or dark matter scalar will not be regarded as an 
independent field, but is induced via (|2"3"]) by the non-Einsteinian pieces of the 
general Lagrangian L — L(R). Solving for the potential via the method of Hclmholz, 
we obtain 

~Rn 



U{tt) = H{R)ltt 2 



2 k 



L 



/Q 2 



(25) 



If we identify the conformal factor with the field momentum via 57 = 2koIL / 'dR, 
the bracket in (I25p can be regarded 35 as a Legendre transformation L — > H(R) = 
RdL I dR — L from the original Lagrangian ^ to the general nonlinear curvature 
scalar Lagrangian L — L(R). Then, the parametric reconstruction 



R = 2«:exp (VW 3< /> 



dU 



2U{cj>) + v / 3/2k — 



(26) 
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Fig. 1. The qualitative form of a Whitney surface with local extrema in the case of a nonlinear 
curvature Lagrangian L(R) depending on </> as control parameter. 



and 



L = cxp (2^2k/3<?(> 



dU 



(27) 



of the higher-order effective Lagrangian L(R) from the self-interacting inflaton 
potential U(4>) arises. 4 Here the scalar field plays merely the role of a control pa- 
rameter. The form of the Whitney surface with its local valleys and mountains is 
qualitatively shown in Fig. [T] 



5. Bifurcations with effective 'dark energy' 

In order to model self-interacting dark matter, let us consider the potential 

U = m 2 \4>\ 2 {l-xW) , (28) 

where m is the mass of an ultra-light scalar and \ the coupling constant of the 
nonlinear self-interaction. It provides us with a solvable non-topological soliton type 
model of dark matter halos 43 ' 44 even with toroidal substructures 45 and a reasonable 
approximation of the rotation curves of dark matter dominated galaxies. Moreover, 
the predicted scaling relation 11 fits almost ideally astronomical observations. In 
Ref. 54 we predicted the effects of such scalar field halos for microlensing. 
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5.1. Free field 

In the string landscape, 25 a triple unification of inflation, DM and DE suggests 
itself based on a simple quadratic potential ([28)) with \ = corresponding to a free 
massive field. Then, our scalar field toy model likewise incorporates 'dark energy' 
in a rather novel way: The exact parametric solution of the equivalent nonlinear 
Lagrangian L(R) for the 'free' field reads 

R = 6m 2 xe x (l + x), (29) 
L = '-^—xe 2x (2 + x), (30) 
where x := lnf2 under the reality condition £1 > 0. 

0.2 
. 1 



-0 . 1 
J-0 . 2 
-0 . 3 
-0 . 4 
-0 . 5 

-1 -0.5 0.5 1 1.5 2 
R 

Fig. 2. The swallow tail behavior of the Lagrangian L(R) for a massive 'free' field (x = and 
K = m = 1). 

The dependence of L(R) given in Fig. [2] is rather surprising and represents the 
bifurcation set of the swallow tail catastrophe 2 ' 21 ' 23 associated with some higher 
dimensional grand manifolds. The scalar field 4>, its mass m as well as x pl av the 
role of control parameters. According to the theory of singularities (more widely 
known as catastrophe theory, cf. Arnol'd 2 ), this bifurcation set indicates that the 
Lagrangian manifolds are associated with two local "minima" and one "maximum" 
(and saddle points at the meeting points of the grand manifold) . Each of the "min- 
ima" merges with the "maximum" at the cuspoidal points A and B and then dis- 
appears. The minima (the semi-infinite segments A and B) are characterized by a 
positive second derivative of the Hamiltonian 7i(f2) with respect to the momentum 
J7, i.e. d 2 TL/dVl 2 > 0. They correspond to an effective Lagrangian with vanishing 
cosmological constant. 

For the "maximum", i.e. the segment AB, this derivative is negative and the ef- 
fective Lagrangian has a modified gravitational constant and a positive cosmological 
constant, 

fveff = en , A off = 3m 2 /(2e) , (31) 
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respectively. 




Fig. 3. Whitney surface for the nonlinear curvature Lagrangian L(R) corresponding to the model 
(1 28 I t with x = 0-3 and m = k = 1 depending on <f> as control parameter. (</> in units of ^/3/(2re) .) 



5.2. Self-interacting scalar 

In the nonlinear case, the exact parametric solution L(R) reads 



R = 6m 2 xe x 



1 



L = 



3m 2 
~2k 



" An?* ' 

_27x 4 
X 2n 2X 



9x , 
4k 2 

4k 2 



(32) 
(33) 



The resulting extremal curve of the Whitney surface is drawn in Fig. [31 and its 
projections L = L(R), R = R(cf>), and L = L(<f>) in Fig. [4] 

For x 7^ 0, the bifurcation diagram in Fig. [5] is of higher rank than for the free 
field. Near the center, we find the butterfly part of the catastrophe. But, for each 
non-vanishing x-, there exists a further cusp far away from the center, for a large 
negative R value and L close to zero. From this cusp, the curve finally returns to 
the origin, for L values very close to zero so that it cannot be seen in Fig. [5] or O 
In an enlargement of the additional cusp in Fig. [71 one can see that L{R) actually 
becomes negative and forms a cusp. 

The effective strength of gravity and the value of cosmological constant 



A 



off 



3m^ 



x + 3 ' 



(34) 
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Fig. 4. Projections of the Whitney surface in Fig.[3]for the nonlinear curvature Lagrangian L(R) 
corresponding to a nonlinear scalar ([28} with \ = 0.3 and m = K = 1. 




5 10 15 20 25 30 
R 



Fig. 5. The wigwam singularity for the nonlinear curvature Lagrangian L(R) for a NTS model 
with x = 0, 0.1 (long dashes), 0.2 (short dashes), 0.3 (dashed) and k = m = 1; here only the 
butterfly part near the center is shown. 

respectively, depend now on both, the mass m and, via x as a solution of R — 
0, on the coupling constant x °f the scalar field, and need to be constrained by 
cosmological data. 

In quintessence models, 65 e.g., the crossover scale of the scalar field is |0 C | = 
exp(— l/a)A/|[/m, where a = 1/138 is about Sommerfeld's fine structure constant 
and Mpi = the reduced Planck mass. Then the tiny observed cosmological 

constant of the present epoch 

A DE = m 2 (j)l = cxp(-2/a)Af^ ~ (l(n 3 eV) 4 (35) 

is roughly reproduced for small x- 

It is important to stress that within the range limited by cuspoidal points, all 
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-20 -10 10 20 30 

R 



Fig. 6. For the NTS model with \ = 0.3 and K = m = 1, the wigwam catastrophe of the 
Lagrangian L(R) can be decomposed into a butterfly plus additional elementary catastrophe. 
Here, the butterfly catastrophe arises near the origin; at about (R, L) = (—250, —0.015) there is 
a further cusp, cf. Fig. [7] the curve finally returns to the origin which cannot be seen due to our 
scale. 



0.04 
0.02 

—I ^ 

-0.02 
-0.04 

-240 -220 -200 -180 -160 
R 

Fig. 7. Continuation of Fig. [6] Enlarging the additional cusp for the nontopological soliton field 
with x = 0-3 and K = m = 1. 

three states may coexist with each other. Thus our bifurcation set, cf. Ref. 55 for 
a generalization with x ^ 0, indicates that each local patch of the Universe may 
have a different strength and cosmological constant ('dark energy') controlled by 
the mass m of the scalar field, but the effective Lagrangian has approximately the 
same Einstcinian form (fT]). In "maximal" domains, inflation may be still going on, 
thus realizing prospective ideas of Linde. 26 

5.3. Natural inflation from the axion? 

In QCD, after integrating out the fermion fields, its generating functional including 
a topological Pontrjagin term for the Yang-Mills gauge fields induces an effective 
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axion potential 

U = A 4 qcb [1 - cos(0// )] . (36) 

This potential displays a periodicity with a period of 2nfy, has a minimum at 
(f> = 0, as required, and leads to the induced axion mass of = Aq CD //0. A 
quintaxion may also be induced by spacetime torsion. 39 

Within natural inflation, 51 such a potential has been proposed for an axion 
coupling constant close to the Planck scale. For simplicity let us assume in the 
following that ~ \/3/ (2k). Then, following the general prescription (l^6]) - ([2T)l . 
we find the reparametrized solution 

R = 2k exp (x) A^ CD [2-2 cos (a:) + sin (x)] , (37) 
L = exp (2x) Aq CD [1 - cos (x) + sin (x)] . (38) 

Fig.[8]exhibits a swallow tail behavior similar as for a NTS potential with x = Q. 

5000 


_, -5000 
-10000 
-15000 

-1 -0.5 0.5 1 1.5 2 
R 

Fig. 8. The axion gives rise again to a swallow tail catastrophe in the Lagrangian L(R), with the 
choice of Aqcd = 0.1, K = 1, and f^ = 1. 




6. Domain structure of dark matter and induced dark energy 

Baryonic particles (such as protons and neutrons) account for visible matter in the 
Universe, that is only a small fraction of the observed total matter. The major part 
of it is a mysterious DM, which only interacts gravitationally. On the other hand 
according to our model, the Universe splits into almost Einsteinian domains with 
depending on the physical scale and eras - a different gravitational strength and 
cosmological constant, the latter being most likely a representation of DE. A natural 
question is how these domains arise? To answer this, we have to look into the initial 
inflationary stage of the Universe, long before decoupling of radiation and matter, 
even before grand unified phase transitions where only an initial pre-field existed. 
At the initial stage of the quantum era, the fundamental pre-field may have existed 
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in a form of a scalar or vector gauge fields. Without lack of generality we can limit 
ourselves to a scalar field <j>. 

During this quantum era, the scalar field was strongly fluctuating. Due to the 
rapid inflation, the size of the Universe as well as the size of each individual fluc- 
tuation have increased enormously while the amplitude remains almost the same. 
During this fast spacelike expansion, just after the stretching, these space fluctu- 
ations of <j) were frozen while loosing their quantum nature due to the increasing 
size. The decoupling of radiation and visible matter has imprinted fluctuations on 
the electromagnetic spectrum which we observe now in the anisotropy microwave 
background radiation. However, we cannot see such ripples of the visible matter 
after the decoupling. Due to the gravitational interaction the process of large scale 
formation has started just after the decoupling. And the state of nearly homoge- 
neously distributed visible matter existing after the decoupling has eventually been 
transformed into clusters of galaxies and star structures which we see nowadays. Ob- 
viously, the decoupling of the dark part of matter from radiation and visible matter 
may have happen even earlier, when the Universe was in the form of a quark-gluon 
plasma. However, after decoupling the dark part of matter must be subjected to 
the same gravitational forces and to a similar large scale formation. Since its total 
mass is significantly larger than the visible mass, the formation of large structures 
of DM may have a weak influence on the visible matter. 

Dark matter will have an opposite effect by dictating where and how the visible 
matter should be distributed, probably localized around clusters of DM. In the case 
of scalar field DM, its peculiarities such as self-interaction may also have contributed 
to the DE distribution. DM cannot be seen directly by traditional observations but 
can be inferred from more sophisticated gravitational lcnsing, 28 revealing that DM 
exists in a form of a loose network of filaments, growing over time, which intersect in 
massive structures at the locations of clusters of galaxies. This finding is consistent 
with the conventional theory of large scale structure formation, where a smooth 
distribution of DM collapses first into filaments and then into clusters, forming a 
cosmic scaffold, 26 then accumulating visible matter and later newly born stars. 

The primary candidate for DM is a scalar field. If a scalar field exists, then it 
may have different amplitudes or different mass densities in different space regions. 
Such a distribution of <f> has been created after the inflation of the Universe. It 
is associated with different amplitudes of the frozen fluctuations. The process of 
large scale structure formation is rather similar in each of these regions. Our re- 
sults indicate that in each of these different regions the gravitational constant may 
vary, 55 depending on the amplitude of the scalar field of the given frozen space 
fluctuation. In other words, in each of these regions defined by the original frozen 
fluctuations there will arise different Newtonian potentials. This is consistent with 
recent studies 9 showing that, during expansion, a vector field having two different 
initial amplitudes will give rise to different Newtonian potentials. This difference 
in turn drives self-consistently enhanced growth of the density perturbations which 
enables tiny perturbations to grow into the large structures we see today. 
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6.1. Cosmic domains 

Our studies indicate that these different Newtonian potentials correspond to differ- 
ent effective gravitational constants in various regions of the Universe. If this is true, 
the data of recent observations 28 should be reinvestigated. In regions with larger 
gravitational constants one may see that the gravitational lensing will be stronger 
and vice versa. We regard areas with different gravitational potentials as cosmic do- 
mains characterized by their own gravitational and cosmological constants. There 
occurs a coexistence of these domains (which primordially may be of topological 
origin 18,19 ). On the boundaries of the domains, changes in the gravitational and 
cosmological constants are expected to be abrupt. In spite of the random or fluctu- 
ational origin of these domains, the values of these gravitational and cosmological 
constants take only a few universal numbers, corresponding to only a finite number 
of types of domains. 

For small curvature, each of these states can be described approximately by the 
same effective Lagrangian (TTJ) with different gravitational constant n e g and cosmo- 
logical constant A e ff and emerge as a fixed point of the conformal transformation 
from one side and universal classes of the smooth differential mappings from the 
other side. 21 These spaces are approximately Einstein spaces, but of different grav- 
itational strength as well as with different cosmological constants. The distribution 
of such cosmological constants in the Universe corresponds to DE induced by the 
scalar field distribution. On a very large scale, (f> is homogeneously distributed and 
therefore we expect that DE will be distributed homogeneously, depending on the 
details of inflationary cosmology. 

6.2. Domains arising from free fields 

If DM and DE were associated with a massive free scalar field gravitationally coupled 
or/and associated with axions (see, Sec. 5.3) there will arise three types of domains: 

1) The first type of domains, I, is described by Einstein's theory with conven- 
tional gravitational constant and vanishing cosmological constant. 

2) The second type of domains, II, is described by an Einstcinian model having a 
very large gravitational constant and a vanishing cosmological constant. Such type 
of domains may be gravitationally unstable, leading eventually to their contraction 
and ultimately to a gravitational collapse and phenomena resembling supernova ex- 
plosions. For very large masses, we may expect that these explosion will be stronger 
than standard supernovae. Our approach may explain the fact why some supernovae 
progenitors seem to have exceeded 15 the Chandrasekhar limit. 

3) The third type of domains, III, are described by an Einstcinian model with 
large gravitational and positive cosmological constant A 1 ^ given by (f34|) . Then we 
expect that such domains are still expanding according to effective Friedman equa- 
tions and probably their sizes are also increasing. Such expansion may be ended 
abruptly and such domains will be transformed into domains of type I. 

One may speculate that the largest part of the Universe is occupied by domains 
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of type I, corresponding to spaces with zero cosmological constant. The next largest 
proportion of the Universe is expected to be occupied by domains of type III. Those 
parts of the Universe occupied by domains of type II are decreasing with time. Each 
disappearance of such domain would be accompanied by supernova type explosions. 
On the average, the cosmological constant inferred from recent observations is 

~ Vi + W+Wi m 

where the spacelike volumes V\, V\\ and V\\\ are associated with domains of type I, 
II, and III, respectively. 

In summary, our bifurcation set (see Fig.[T]and Fig. [2j) indicates that the Universe 
is locally described by Einstein's GR and effectively split into domains. The splitting 
originated from primordial quantum fluctuations which were frozen during inflation. 
Each domain or local patch of the Universe may have different gravitational strength 
and may have zero or negative pressure associated with DE. In domains of type II 
associated with a positive cosmological constant the inflation may be still going 
on. It will be stopped exactly at the bifurcation point A and this domain will be 
transformed into a domain of type I, as in Fig. [5J For a massive free scalar field and 
for axions the bifurcation set is the swallow tail catastrophe. There are two cusps 
at the points A and B which are associated with the highest singularities of the 
differential mappings. They are also related to the domain boundaries. In each of 
these cuspoidal points A and B of the bifurcation diagram, the minimum of some 
grand manifold merges with the maximum. 22 Each local minimum (the segment A 
and the semi-infinite segment B) is associated with two types of domains, I and 
II, respectively. Each local maximum (saddle) of the grand manifold is associated 
with an expanding domain of type III (the segment AB). The effective Einsteinian 
gravity for this domain is stronger than that in the domain of type I, depending on 
the mass of the scalar field. Domains with a positive effective cosmological constant 
are still in an expanding inflationary phase. 

The 'strong' gravity state, the domain II (the semi-infinite segment B) with a 
negative or vanishing cosmological constant corresponds to deflation. Since the Uni- 
verse is split into dynamical regions associated with different gravity, their bound- 
aries change with time and present some kind of cosmological strings or membranes. 
Changes of the boundary are associated with some sort of local phase transition re- 
minding us of phase separations like those arising in first order phase transitions. At 
some instant, all these domains are in quasi-equilibrium with each other although 
their boundaries changes with time as the membranes are moving. At very long 
time scales (much larger than the Planck time) when inflation in the domains of 
type III will eventually be ceased, only stable phases, like the domains of type I 
and some small concentration of domains of type II will remain. These domains 
would produce explosions resembling supernovae with masses well beyond of the 
Chandrasekhar limit, similar to those observed. 15 
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6.3. Domains arising from self-interacting fields 

In case of self-interacting scalar fields with x ^ 0, the classification will be richer 
(see Figures |H [5] and [6]) , with five types of domains. The bifurcation diagram as- 
sociated with such a wigwam catastrophe or ^5 singularity according to Arnold's 
classification 2 consists of four cusps or four A 2 singularities, with a dramatic change 
in the gravitational and cosmological constants. 

There occurs a general decoupling of the wigwam catastrophe A5 into a butterfly 
catastrophe A±, and the elementary cusp A2, depending on the values of R: If the 
butterfly catastrophe arises near the origin or at small values of R, the additional 
cuspoidal point is associated with a very large value of R ~ —250. For some fixed 
values of the coupling constant x = 0.1, 0.2, 0.3, the butterfly part of bifurcation 
diagram is presented in Figs. [5] In comparison with the free massive field, the in- 
clusion of the self-interaction associated with the nonzero parameter \ gives rise to 
two extra cusps leading to the appearance of two extra branches in the bifurcation 
diagram, again associated with effective Einsteinian spaces. Although these almost 
linear branches do exist at positive R, their major part is related to negative val- 
ues of R. For one of these branches the cosmological constant vanishes while for 
the other one the cosmological constant decreases when the value of x increases, 
and the part of the branch associated with R > decreases. In general, the self- 
interaction induced by the field 4> increases the strength of gravity and decreases 
the cosmological constant. 

Let us describe in more detail the transitions in the bifurcations: Starting from 
the first branch, denoted as I, through the first cuspoidal point we arrive at the 
second branch II which is very similar to the case of a free massive scalar field, 
although its length decreases significantly when x increases and the cosmological 
constant vanishes at any value of x- This branch corresponds to the limit <fi — ► 
of small scalar fields, corresponding to the conformal factor = 1. Therefore, on 
this branch we recover the linear Hubert-Einstein Lagrangian Lhe = R/2k, as ex- 
pected. With decreasing fi, we will come to the next cuspoidal point, where the new 
branch III has a positive cosmological constant and gravity becomes stronger. With 
increasing value of x the length of this branch increases. Thus, there will be more 
domains having positive cosmological constant and a gravitational strength higher 
than the conventional. By decreasing Q until reaching the next cuspoidal point of 
the bifurcation diagram, we arrive at the branch IV. The spacetime associated with 
this branch has a very strong gravity and negative cosmological constant. Similar 
to the free massive case, this branch probably is gravitationally unstable. Finally, 
for very small f2, we will arrive at the branch V with R < 0, where gravity is the 
strongest but without cosmological constant. Similar to the previous branch IV, the 
domain V is unstable, as well. Thus in a Universe filled with self-interacting scalar 
fields, we would expect that a rather strong phase separation into further domains 
will arise. 
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7. Discussion 

For dilute DM the Lagrangian has the standard Hilbert-Einstein form Lhe = R/2k. 
For dense DM, i.e. large values of the DM scalar, the resulting effective Lagrangian 
is Lhe = R/2k c s, only the slope is less steep, i.e. the effective gravitational coupling 
n c ff > k is, by "renormalization" , 3 ' 65 larger than Newton's. This resembles some 
aspects of MOND, but with the advantage that we still work in the standard general 
relativistic framework. Depending on its sign, the cosmological constant A c g in a 
bifurcation approximated for small by (JTJ) can also model DE or an accelerating 
phase of the present epoch of the Universe (like 'anti-gravity'). There also may 
arise "gravitational screening" for K e s smaller than Newton's. MACHOs can now 
be described by some specific nonlinear modifications of the Hilbert-Einstein action. 

In the light of our finding we have to comment on recent observations, 29 which 
constrain the density of MACHO type objects in the Universe by measuring the 
brightness of high redshift type la supernovae relative to low redshift samples. These 
data favor DM made of microscopic particles (such as WIMPs) or scalar and vector 
gauge fields over MACHOs with masses between 1 -2 and 10 10 solar masses. This 
provides another evidence that our approach correctly describes the frozen space 
fluctuations which lead to a formation of cosmic domains with a different density 
of DM and DE in the Universe. 
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